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Abstract 

The nuclear incompressibility is deduced from measurements of the Isoscalar Giant Monopole 
Resonance (ISGMR) in medium-heavy nuclei, and the resulting value turns out to be model de- 
pendent. Since the considered nuclei have neutron excess, it has been suggested that the model 
dependence is due to the different behaviour of the symmetry energy in different models. To clarify 
this issue, we make a systematic and careful analysis based on new Skyrme forces which span a 
wide range of values for K^o , for the value of the symmetry energy at saturation and for its density 
dependence. By calculating, in a fully self-consistent fashion, the ISGMR centroid energy in ^'^^Pb 
we reach, for the first time within the non-rclativistic framework, three important conclusions: 
(i) the monopole energy, and consequently the deduced value of -f^oo, depend on a well defined 
parameter related to the shape of the symmetry energy curve and called Kgym', (ii) Skyrme forces 
of the type of SLy4 predict Koo around 230 MeV, in agreement with the Gogny force (previous 
estimates using Skyrme interactions having been plagued by lack of full self-consistency) ; (iii) it is 
possible to build forces which predict K^c around 250 MeV, although part of this increase is due 
to our poor knowledge of the density dependence and effective mass. 
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I. INTRODUCTION 



The question about the proper value of the nuclear incompressibility is still open. 
The model dependence of this quantity amounts to a difference of the order of ~ 10-20% 
among the values obtained within different theoretical models. There is a renewed interest in 
this issue, motivated both by the improved quality of the recent experimental measurements 
of the Isoscalar Giant Monopole Resonance (ISGMR), and by the progress of Relativistic 
Mean Field models (RMF) which are to be confronted with more traditional non-relativistic 
models based on Skyrme or Gogny effective forces. 

Skyrme energy functionals have been widely used in nuclear structure calculations over 
three decades. The first Skyrme effective forces were built in the pioneering work of Vau- 
therin and Brink [1], by fitting their parameters to nuclear matter properties (the saturation 
point) and to selected observables (binding energies and charge radii) of closed-shell nu- 
clei calculated in the Hartree-Fock (HF) approximation. Later, many improvements of the 
Skyrme energy functional have been devised. These have been possible also because the 
mean field approach was extended to the time-dependent case (time-dependent HF, TDHF) 
and to its small amplitude limit (Random Phase Approximation, RPA). Within this scheme, 
it is possible to calculate the collective nuclear excitations and to explore the correlations 
between their properties and the force parameters, or some physically meaningful combina- 
tions of them. The relation between the ISGMR and the nuclear incompressibility is one of 
such relations. 

The introduction of reliable RMF effective Lagrangians is more recent. However, the 
progress in this field has been quite fast [2] , and we can nowadays discuss the properties of 
the RMF parametrizations on the same footing as the Skyrme [3] and Gogny [4] functionals. 

The nuclear incompressibility K^o is related to the curvature of the energy per particle 
E/A in symmetric nuclear matter around the minimum qq, i.e., at the saturation point: 

Ko. = 9gl-^j\so- (1) 

The interest of determining the value of stems also from its impact on the physics of 
supernovae and neutron stars. 

Until a year ago, the status of the nuclear incompressibility problem could be summarized 
as follows. From calculations based on Skyrme functionals different groups have extracted 
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values of K^o of the order of 210-220 MeV. Using the Gogny functionals, a value of K^o 
around 230 MeV was obtained. Finally, the relativistic calculations predicted values in the 
range 250-270 MeV. All these results made use of the measured value of the ISGMR, e.g. 
in ^°^Pb, as explained below. 

This situation for the nuclear incompressibility "puzzle" has been reviewed in Ref. [5]. 
There, it was shown that the accuracy of the ISGMR data obtained at Texas A&M [6] 
allows for extracting K^x, with an experimental error of no more than ± 12 MeV. Moreover, 
a rather important new conclusion was reached. The previous works based on the Skyrme 
forces consisted of not fully self-consistent HF plus RPA calculations in which the two-body 
residual Coulomb and spin-orbit forces were neglected. These terms are rather small since 
they affect the monopole energy in ^osp]-, only 4-5%, but this produces a change of 8-10% 
in the extracted value of the nuclear incompressibility. By considering this effect, the new 
value of from Skyrme functionals turns out to be 235 MeV, in very good agreement with 
that extracted from the Gogny calculations. Consequently, there is no discrepancy between 
the results of the different non-relativistic calculations. On the other hand, the gap with the 
relativistic results remains significant. 

The most recent attempts in the literature [7, 8, 9, 10] to attack this problem are focused 
on the possible relation(s) between the monopole energy in systems with a neutron excess, 
hke ^°®Pb, and the density dependence of the symmetry energy S{g). In fact, one of the 
clear differences between the Skyrme and RMF functionals concerns the behaviour of the 
symmetry energy around the saturation point ^o- The Skyrme energy functionals are char- 
acterized by smaller values of the symmetry energy at saturation, and of the corresponding 
slope, as compared with the RMF functionals. In this sense it may be said that the RMF 
functionals are "stiff" compared with the "soft" non-relativistic ones. 

In Ref. [7] some effective Lagrangians whose symmetry energy has different density depen- 
dences are built. This is easy to achieve since, by adjusting the p-meson coupling constant, 
one can at the same time soften the symmetry energy S{g) and lower its value at the sat- 
uration point, J = S{go)- It is thus found that the extracted values of K^o indeed differ 
and can even become close to the Skyrme force values if J is around 28 MeV. However, in 
Ref. [7] no systematic treatment of finite nuclei is attempted. In Ref. [8] it is pointed out 
that RMF parametrizations with J lower than 32 MeV cannot describe satisfactorily the 
N Z nuclei. The authors conclude from their calculations that the lower limit for the 
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RMF value of K^o is around 250 MeV. In Ref . [9] , using a markedly improved version of the 
model of Ref. [7] , this lower limit is confirmed since K^o results to be 248 MeV. 

While in the relativistic framework it seems impossible to push the value of below 
this lower limit, the recent results of Ref. [10] suggest that one can build at least one Skyrme- 
type interaction having K^=2bb MeV and reproducing the correct ISGMR energy in ^°^Pb. 
This is at variance with all other non- relativistic calculations quoted in Ref. [5]. Moreover, 
the origin of the result of Ref. [10] is unclear. Does it correspond to an effective force with 
a symmetry energy which is as stiff as that associated with the RMF Lagrangians ? Or is 
some other parameter playing a role ? 

It seems necessary to make a systematic analysis of what is the upper limit for K^o 
within the non-relativistic framework. The value of 230 MeV is extracted from a subset of 
the existing Skyrme parametrizations. Our main goal is to answer the question whether it is 
possible to build new interactions in such a way that becomes closer to the RMF value, 
as well as to study what are the crucial quantities which control this variation of i^oo- In 
particular we would hke to understand the (so far, singular) result of Ref. [10]. 

The structure of the paper is the following. In Sec. II we review how the nuclear incom- 
prcssibility is extracted from the microscopic calculations using the ISGMR experimental 
data, and what are the plausible quantitative arguments in favour of the idea that the den- 
sity dependence of the symmetry energy plays a role. In Sec. Ill we describe our fitting of 
new Skyrme interactions suitable for the microscopic ISGMR calculations, and in Sec. IV 
we describe the results obtained and the implications for the value of i^oo- In Sec. V we 
present our conclusions. 



II. DEDUCING THE NUCLEAR INCOMPRESSIBILITY FROM MONOPOLE 
DATA AND THE ROLE OF THE SYMMETRY ENERGY 

In all the discussions of the relationship between the nuclear matter incompressibility and 
the ISGMR in finite nuclei, the starting point is the definition given by J.R Blaizot [11] of 
the finite nucleus incompressibility as 

-^A= ^2 ' 
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where m is the nucleon mass and (r^)o is the ground-state mean square radius. This expres- 
sion has a well-defined meaning in medium-heavy nuclei, where the ISGMR is associated 
to a single peak at the energy £^isgmr ~ SOA~^^^. In light nuclei the monopole strength 
is very much fragmented and many states show up, whose microscopic structure does not 
correspond to the simple picture of the radial "breathing mode" according to theoretical cal- 
culations (see, for example, Ref. [12]). In the case of the nuclei studied in [6], the existence 
of a single, collective monopole state is quite evident from the measured cross sections. In 
particular, in the case of ^°^Pb which is the object of our present study, the experimental 



peak energy and the centroid energies Eq and E^i (defined respectively as ^ and 



where is the k-th moment of the strength function) essentially coincide, leaving out any 
ambiguity about the correct value of £^isgmr to be used for determining the experimental 
value of i^A- 

However, finding a theoretical relation between Ka and is less simple. In Ref. [11], 
the generic expression of the energy functional associated with Skyrme-HF has been written 
in the case of a finite spherical system. At variance with that of infinite matter, the density 
is not uniform and cannot be reduced to a simple number. Therefore, to minimize the energy 
functional and find its second derivative around the minimum one has to resort to various 
simplifying hypotheses. The main one is the use of the so-called scaling model, in which 
a simple shape of the ground-state density is assumed and its changes are associated 
to a single parameter A, i.e., they are of the type go{r) Qx{r) — -^Qoij)- In this way, 
the expression for the finite system incompressibility can be found. By isolating the terms 
corresponding to the volume, surface, symmetry and Coulomb contributions, the result can 
be written as 



We remind here that, in the past many authors have used the formula (3) as an ansatz 
and have tried to obtain the parameters of the r.h.s. from a numerical fit, using as input 
the experimental values of K\ in different nuclei. This procedure is not stable and leads to 
ill-defined values of the parameters [13], so that it is nowadays abandoned. 

Instead, the microscopic method to deduce K^o rehes on the fact that RPA calculations 
of the ISGMR can be performed by using functionals characterized by different values of 
Koo- If the calculations done with a given functional reproduce the experimental ISGMR 




Coul 



A4/3' 



(3) 



where S={N-Z)/A (cf. Sec. 6.2 of Ref. [11]). 
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energy, the associated value of K^o should be chosen as the best one. Let us examine this 
in more detail. 

Mainly one nucleus has been used so far, that is, ^°^Pb. In the first work in which the 

microscopic procedure has been applied [14], the RPA values for Kp^ obtained from the RPA 
centroid energies E_i have been plotted versus the of the force used. Then, an empirical 
linear fit of the results was performed, namely 

i^A = aK^ + h. (4) 

This relation allows to extract the best value for Kao by inserting the experimental K^. In 
[14] the explicit form of (4) in the case of ^ospb is Ka = 0.647^00 - 3.5 [MeV]. The second 
term of the r.h.s. is much smaller than the first term. Consequently, even if in principle the 
last formula together with (2) would lead to £'isgmr = 1.16-\/0. 64/^00 ~ 3-5, this equation 
can be approximated by 0.928-\/^oo (neglecting the second term under the square root). 
This explains why in many of the works quoted in [5] a successful interpolation of the type 

EjsGMR^a'^fKZ + b' (5) 

was done: in practice, Eqs. (4) and (5) are equivalent. It is from either of these relations, 
using the experimental ISGMR energy in ^°^Pb which is 14.17 ± 0.28 MeV, that the values 
for Kfx, mentioned in Sec. I were obtained. 

The uncertainity of the value of which is deduced is 

SKqo _ ^ 6EisGMR 

The experimental error on the monopole energy, plus a theoretical error of the same order 
(see [5]), produce a global error bar of ± 12 MeV on K^o- 

One may argue why the linear relations just introduced are valid. So far, Eq. (3) has 
not played in fact any explicit role in the deduction of K^o- However, this expression can be 
taken as a rather useful guideline. Given a microscopic functional, the different terms K^urf, 
Ksym and Kcoui (in addition to i^'oo) entering this formula can be calculated as described 
shortly below. The resulting value of differs from the microscopic outcome of RPA, as 
a rule, by about 5%. Therefore, we make in the rest of this Section a detailed analysis of 
the role of the different terms in (3). If, for a family of functionals, Kp^ depends linearly 
on as written in (4), it means that the other terms do not vary significantly. This is 
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what happens for a large subset of the Skyrme and Gogny parametrizations, as it is evident 
from Fig. 6 of Ref. [14]. However, the role of the surface, symmetry and Coulomb terms 
should be critically re-examined if new functionals, including relativistic ones, enter in the 
discussion. 

The expression for these terms have been given in Ref. [11]. K^urf cannot be calculated 
analytically, but numerical estimates are possible within both the quantal and semiclassical 
scheme. We refer the reader to Ref. [15] for an example of a quantal derivation (which is 
a scaled Hartree-Fock calculation of semi-infinite nuclear matter). The most recent semi- 
classical, i.e., extended Thomas- Fermi (ETF) calculations, have been performed both in the 
non-relativistic and in the relativistic scheme and have shown that the quantity Kg^rf is well 
approximated by cKoo with c ~ -1 (however, it should be noted that c tends to grow with 
Koo) [16]. We have checked that this approximation is valid in the case of all the forces used 
in this work: we have seen that, e.g., c = —1.03 if Xoo^SSO MeV and c = —1.07 if ^00=250 
MeV. 

In order to study Kgym, we first give some necessary definitions of the symmetry energy 
and of the parameters related to its density dependence. We define the symmetry energy 
by writing the total energy density S as the sum of an isoscalar part So{g) which depends 
only on the total density g = gn + Qp, and an isovector part, 

£ig,g.)=£o{g) + gSig){^f, (6) 

Q 

where g_ = g^ — gp. We remind in this context that in a homogeneous system, E/A = S/g. 
The symmetry energy S{g) can be expanded up to second order around go, 

S{g) = S{go) + S'{go){g - go) + ls"{go){g - gof. (7) 

The value of the symmetry energy at saturation S{go) is often denoted as J and we are 
following use the same notation in this paper. Other notations, like o-^ or 04, are also 
employed in the literature. The first and second derivatives of S{g) at the saturation point 
have been written many times in terms of the so-called parameters L and Kgy^ (see, e.g., 
Ref. [17]), as S'{go) = L/3go and S"{gQ) = Kgym/Qgo- It is quite unfortunate that the symbol 
Kgym has been used in the literature with such different meanings, either in connection with 
S"{go) or in Eq. (3). Here, we will always use Ksym to mean the symmetry term of K\ in 
Eq.(3). 
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The expression of Kgym is 

Ksym = 9glS"{go) + 9goS'{go) - ^^^^^^U, (8) 

and from this expression it is evident that this parameter contains some relevant information 
about the density dependence of the symmetry energy. 

The values of J rule, larger in the case of the RMF functionals than for the 

Skyrme ones. A larger value of J is correlated with a larger value of 5"(^o), which is usually 
a positive quantity although it may sometimes become negative (cf. [18] and Fig. 4 of [8], 
as well as Tables 1 and 11 in this paper). In the next Section we show that a larger J is also 
correlated with a more negative value of Kgym, at least for the forces we have studied. The 
explanation which is given for the correlation between J and S'{go) is that the fits to finite 
nuclei observables constrain the symmetry energy at some average density (g) lower than 
^0 (see, e.g., Ref. [19] and references therein). In the case of one set of forces introduced 
in this paper (see Sec. Ill), this typical behaviour of the symmetry energy is shown in Fig. 
1. In a narrow region around ^=0.10 fm~^ (± 0.001 fm~^) all curves cross one another at 
a value S{g)=25 ± 1 MeV. When the symmetry energy at saturation is larger, the slope is 
also larger. The other sets of forces show qualitatively the same trend. 

The last term of (3) is the Coulomb contribution which is unlikely to be very much model 
dependent. It is written as 

where ro is the average inter-particle spacing. 

In summary, if we want to compare two models, say I and II (they could be for instance a 
non-relativistic and a RMF functional, respectively), we will write, by using Kgurf — ci^oo 

^ Kg'\l + CA-V3) + ^(11) 52 + 41;)^^. (10) 

We have already mentioned that Kgym is negative, and the same is true for Kgurf and Kcoui- 
All can be viewed as corrections to the leading term K^o- It is clear that a more negative 
value of K surf or Kgym leads to extracting from the experimental a larger value of K^. 
We will develop this argument in Sec. IV. 
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III. CONSTRUCTION OF NEW SKYRME PARAMETER SETS 



The different forces used in this study have been built using a procedure which is quite 
similar to the one discussed in Ref . [20] . The starting point is the standard form of a Skyrme 
interaction as given in Eq. (2.1) of [20]. 

In the case of the first set of forces that we have constructed, the density-dependent term 
has = g^^^. The spin-gradient terms occuring in the Skyrme functional are neglected and 
the Coulomb exchange term is included within the Slater approximation. The center-of-mass 
motion is taken into account with the usual correction in the kinetic term, which means 
that only the one-body part of the center-of-mass energy is subtracted before variation. 

The parameters of the forces have been determined by minimizing a built on: 

1. the infinite nuclear matter properties ^o, E/A{qq) (while Xoo, J and the enhancement 
factor K of the Thomas- Reiche-Kuhn sum rule are kept constant); 

2. the following finite nuclei properties: binding energies and charge radii of ^°'^^Ca, ^^Ni 
and ^°^Pb together with the binding energy of ^^^Sn; 

3. the spin-orbit splitting of the neutron 3p shell in ^°^Pb; 

4. the surface energy, calculated in the ETF approximation and fitted to the value of 
the SkM* force in order to obtain good mean field properties at large deformations 
(especially a good fission barrier of the ^^°Pu nucleus). 

Furthermore, the parameter X2 is fixed to —1 in order to ensure the stability of the fully 
polarized neutron matter in a simple but tractable way [21]. Unlike the case of the SLy4 
force, the equation of state of neutron matter is checked but not fitted in order to have 
a large enough variational space of parameters when the nuclear incompressibility and the 
symmetry energy are varied. The forces which have been built have equal to 230, 240 
and 250 MeV whereas J is varied between 26 and 40 MeV. In Fig. 2 we show the accuracy of 
the present forces in reproducing the ground-state observables (binding energies and charge 
radii). 

Motivated by the comparison with Ref. [10], we have also built another set of forces 
with a similar protocol but with the density dependent term q"' having the same exponent 
q;=0.3563 as the force SK255 introduced in Ref. [10]. The forces of this set have equal 
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to 250, 260 and 270 MeV while J is varied between 28 MeV and 42 MeV. Fig. 3 gives an 
idea of the accuracy of this set, in the same way as for the previous one. 

The nuclear matter properties associated to all the new forces introduced in this paper 
are summarized in Tables I and II. By looking at the values of the effective mass, one can 
recognize the well-known correlation between K^, a and m* jm (see, e.g.. Fig. 2 of Ref. 
[20]). 



IV. RESULTS AND DISCUSSION 



Using these new Skyrme interactions, the ISGMR centroid energies — in ^'^^Pb 

have been calculated in a fully self-consistent manner. The energy-weighted sum rule mi 
is obtained from the well-known double commutator expectation value, while the inverse 
energy-weighted sum rule m_i is extracted by means of a constrained HF (CHF) calculation 
[22] . Adding to the Hamiltonian a term AM, where M is in this case the monopole operator 
YH=irf and A > (to avoid an Hamiltonian without lowest bound), the value of m_i can 
be extracted in two different ways, that is, 

^-^ = -2^ = 2^A^- 

By varying the steps in A and by comparing the outcome of these two different expressions, 
numerical tests concerning the accuracy of m_i can be performed. We have come to the 
conclusion that this quantity can be determined with an accuracy of ± 3% or better. This 
is definitely more reliable than the result of usual RPA calculations made using a basis 
expansion, since the convergence of m_i with the basis size can be quite slow. Moreover, 
as already discussed in Sec. I, Skyrme RPA calculations of the ISGMR performed so far 
lack full self-consistency since part of the residual interaction (the two-body Coulomb and 
two-body spin-orbit terms) are dropped. This has been shown to lead to a systematic error 
in the monopole centroid energies [5] . 

In Figs. 4 and 5 we show the results for the monopole energy obtained with the 
present interactions, as a function of the associated values of and J. Fig. 4 refers to 
the forces with q;=1/6, whereas Fig. 5 is for those with «=0.3563. In the case of the forces 
with q;=0.3563, we have used the same procedure of Ref. [10], that is, we have neglected the 
Coulomb exchange and we have omitted the center-of-mass correction in the HF variation. 
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by subtracting it afterwards from the total energy, in the harmonic oscillator approximation 
{Ecm. — |41A~^/^). We have checked that this lowers the monopole energy by about 150 
keV. The straight hues in the figures are hnear fits of the CHF results corresponding to the 
different symbols, whereas the experimental range for the monopole energy [6] is delimited 
by the horizontal lines. These figures can be compared with Fig. 5 (upper panel) of Rcf. [8]. 
The results for the monopole energy are, as expected, much less sensitive to J than to /Too- 
By varying Koo by 10 MeV, i.e., by about 4%, the monopole energy changes by 0.5 MeV. 
In order to obtain the same change, J should be varied from 26 MeV to 40 MeV which is 
about 50%. The RMF results show quahtatively the same pattern. 

We have to stress that the existence of a definite, yet not strong, dependence on J is 
in agreement with the discussion in Sec. 11, where the role of Kgym as one of the crucial 
parameters governing the monopole energy has been emphasized. It is clear from Figs. 6 
and 7 that the monopole energies do depend on the parameter Ksym, associated with the 
density dependence of the symmetry energy. On the other hand, in the forces we have built, 
there is a strong correlation between Ksym and J, essentially independent of but not of 
a. This correlation is shown in Fig. 8. The interpretation of this plot is the following. The 
modification of the exponent a in the Skyrme functional, allows us to change the value of 
Kgym keeping fixed the values of K^x^ and J. According to the argument developed at the 
end of Sec. II, this should allow to change the value of extracted from the experimental 
ISGMR data. 

By considering only the set with a— 1/6, we confirm the previous result of Ref. [5] that 
i^oo=230-240 MeV is the preferred value for the nuclear incompressibility. This is not fully 
compatible with the RMF result. In fact, extrapolating from Fig. 4, one can see that an 
hypothetic Skyrme parametrization having that associated value of K^o, would reproduce 
the experimental monopole energy only with an unrealistic value of J above 50 MeV. On 
the other hand, the set of forces with q;=0.3563 allows extracting a value of K^o around 250 
MeV, in agreement with the outcome of the RMF calculations. 

We should at this point stress that the more negative values of Kgym which characterize 
the forces with a=0.3563 cannot, alone, explain the extraction of a larger Koo- The forces 
which better reproduce the experimental monopole energy are those with (a, K^, J) given 
either by (1/6, 230 MeV, 28 MeV) or by (0.3563, 250 MeV, 30 MeV). We can apply Eq. 
(10) to the case of these two forces, by inserting the values given in the Tables I and II and 
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by taking into account that c ~ -1 for K^^230 MeV but ~ -1.1 for K^^250 MeV. This 
gives for the two forces, respectively, 

Ka = 154.16 = 230 - 38.82 - 10.23 - 26.79 [MeV] 

and 

Ka = 160.64 = 250 - 46.41 - 15.18 - 27.77 [MeV], 

where the four numbers on the r.h.s. correspond respectively to the volume, surface, sym- 
metry and Coulomb contributions. It is clear that most of the gain of 20 MeV in K^o comes 
from the increase of the surface term (8 MeV), and the more negative Ksym which, mul- 
tiplied by the tiny value of 5^ for ^°^Pb, contributes by 5 MeV. A small contribution of 1 
MeV results from the increase of Kcoui- Finally, Eq. (10) does not consider that, in the 
calculations done by employing the forces with q;=0.3563 the Coulomb exchange and center- 
of-mass corrections are neglected. As mentioned above, this lowers the ISGMR energies by 
about 150 keV and hence Ka by about 5 MeV. This brings the two results for Ka rather 
close to each other. 

V. CONCLUSION 

Until recently, the extraction of the nuclear incompressibility from the monopole data 
was plagued by a marked model dependence: the Skyrme energy functionals seemed to 
point to 210-220 MeV, the Gogny functionals to 235 MeV, and the relativistic functionals 
to 250-270 MeV. It has been shown in Ref. [5] that the result of the Skyrme functionals 
is in fact consistent with that of Gogny, i.e., 235 MeV using the ^osp^ data. The previous 
value of 210-220 MeV was derived using non fully self-consistent calculations neglecting the 
residual Coulomb and spin-orbit interactions. The discrepancy between the non-relativistic 
value of 235 MeV and the relativistic prediction remained, since new relativistic calculations 
confirmed the lower bound of about 250 MeV. 

The work of Agrawal et al. [10] suggests that it is possible to build Skyrme forces which 
fit nuclear ground states and lead to the correct monopole energy, with K^—2hh MeV. Very 
recently, Goriely et al. [23] have also obtained Skyrme parametrizations which fit very well 
the binding energies of a quite large number of nuclei. Some of the parametrizations having 
(a =0.5, m*/m =0.80, J =30 MeV, K^ =276 MeV) are doing satisfactorily well for the 

13 



ISGMR energy in ^°^Pb (-Eisgmr — 14 MeV). In the present work we systematically explore 
what are the conditions which lead to such different results for K^o within the Skyrme 
framework. 

To this aim, we build classes of Skyrme forces which span a wide range of values for 
and for the symmetry energy at saturation J. All these forces reproduce the ground state 
observables with good accuracy We use them to calculate the monopole energy in ^°®Pb, 
defined as E^i — sj^^l- stress again that we can obtain this quantity without any lack 
of self-consistency, and with a numerical error which is not larger than the experimental 
uncertainity 

A first class of forces are built using the SLy4 protocol and have a density dependence 
characterized by the exponent q;=1/6. With these forces, a value of K^^ around 230-240 
MeV is obtained, confirming the previous results of [5]. To obtain the correct monopole 
energy with larger values of would require an unrealistically large value of J, since 
-E'iSGMR increases with K^o and decreases with J. We understand this latter dependence 
as a consequence of the direct relation between the K^o and the density dependence of the 
symmetry energy S{q), and of the unavoidable correlation between S{q) and J. 

To solve the discrepancy with the result of Agrawal et ai, we have built a second class of 
forces which have the density dependence q;=0.3563. Using this class of forces we can arrive 
at Koo between 250 and 260 MeV. Actually, we can reproduce very accurately the results of 
Ref. [10], if we use the same approximations made in that work, namely if we neglect the 
Coulomb exchange and center-of-mass corrections in the HF mean field. This shows that 
the differences between Ref. [10] and our work in the detailed protocol used to determine 
the forces, are unimportant. We have observed that the differences between the results of 
the two classes of Skyrme forces built in the present paper, come both from the surface 
and symmetry contributions, as a consequence of the change in the exponent a, and from 
the neglect of Coulomb exchange and center-of-mass corrections, which affect the monopole 
energy by about 150 keV and therefore, K^o by about 5 MeV. 

In conclusion, within the non-relativistic framework there is not a unique relation between 
the value of associated with an effective force and the monopole energy predicted by 
that force. Bona fide Skyrme forces can either predict 230-240 MeV for or arrive at 250 
MeV if a different density dependence is adopted and if one excludes some terms from the 
energy functional. This latter procedure, although it may mimick the relativistic case, is 
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not conceptually well justified. 
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FIG. 1: Density dependence of the symmetry energy for one of the set of forces (a=l/6 and 
Koo=2A0 MeV). The symmetry energy J is varied from 26 to 40 MeV. 

FIG. 2: Difference between experimental binding energies (left) or experimental charge radii (right) 
with the predictions of forces characterized by a=l/6 and by different values of Koo and J, for 
typical spherical nuclei. Note that the binding energies of ^°'^^Ca, ^^Ni, -"^^^Sn and ^°^Pb, as well 
as the charge radii of ^°'^^Ca, ^^Ni and ^°^Pb are used in the fit of the force parameters. 

FIG. 4: The ^'^^Pb ISGMR centroid energy E^i calculated with the Skyrme parameter sets with 
a=l/6, as a function of J. The different symbols correspond to the values of Kry^ (see inset). 
Lines are numerical fits and are simply intended to guide the eye. The area delimited by the two 
horizontal lines correspond to the experimental value. 



FIG. 5: Same as Fig. 4, for the Skyrme parameter sets with a=0.3563. 



FIG. 6: Same as Fig. 4, with the results plotted as a function of K. 



FIG. 7: Same as Fig. 5, with the results plotted as a function of K, 



FIG. 8: Correlation plot Kgym vs. J. Squares and triangles display respectively the results for the 
forces with a=l/6 and 0.3563. Straight lines are numerical fits. 

FIG. 3: Same as Fig. 2 in the case of the forces with q;=0.3563. 
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TABLE I: Nuclear matter properties calculated with the different Skyrmc parameter sets charac- 
terized by a=l/6 and by different values of Koo and J (these two quantities identify the parameter 
set and are shown in the first column). All quantities are defined in the text. 
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TABLE II: The same as Table I for the different Skyrme parameter sets characterized by q;=0.3563. 
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